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Ferroic domain walls occur naturally as intrinsically nanoscale interfaces separating dif- 
ferent orientations of spontaneous electric polarisation, magnetization, or strain. They can 
present radically different properties from their parent materials, opening an alternative 
route towards unusual functionalities with enormous application potential [TJ [2] . Much re- 
cent work has focused on the complex multiferroic insulator BiFe03 [3j, with discoveries of 
possible magnetism in 109° domain walls [HE], an d conduction in 109°, 180° [6], and sub- 
sequently 71° [7] domain walls. To date, no clear consensus has emerged on the microscopic 
origin of the conduction j6j|8], variously attributed to band gap lowering, or the presence of a 
potential step due to a polarisation discontinuity attracting charged defects. Both scenarios 
relate the conduction to the strongly non-Ising nature of BiFe03 domain walls [9]. Here, 
we demonstrate that 180° domain walls in the much simpler, purely ferroelectric tetragonal 
perovskite Pb(Zr .2Ti .8)O3 (PZT) are also conducting. The observed domain wall conduc- 
tion is clearly differentiate from polarisation switching currents and shows highly nonlinear 
and asymmetric voltage characteristics, strong temperature dependence, and high stability. 

At BiFeOa domain walls, oxygen vacancy defect states inside the band gap clearly modu- 
late the conduction [TO] , but scanning tunneling spectroscopy measurements [H] also show 
effective band gap narrowing at 109° and 71° domain walls. A wide range of conduction 
mechanisms, both bulk- and interface-limited, have been suggested from local probe studies 
[S ED]- To clarify the mechanism of domain wall conduction, measurements in a purely 
ferroelectric, tetragonal perovskite system such as PZT, without the complexities and mul- 
tiferroic nature of BiFeOs, are extremely useful. Interestingly, PZT and PbTiOs show little 
change in band gap values between ferroelectric (3.6 eV) and paraelectric (3.4 eV) phases 
[12], suggesting that band gap narrowing would not play a significant role. Even the "simple" 
180° domain walls in this material can show a complex internal structure, with polarisation 
rotation and flux-closure under specific electric or strain boundary conditions [131 E] • 

Recent transmission electron microscopy measurements beautifully demonstrate such po- 
larisation rotation in a narrow region at the junction of a 180° domain wall and the interface 
of a PZT film and its 1.5-unit-cell-thick SrRu03 substrate [15]. Such features present a 
potential step necessitating screening, which could occur via charged defects, providing 
donor/trap states near the ferroelectric surface and facilitating electron emission/tunneling 
from an adjacent electrode. Alternatively, internal charge redistribution in semiconducting 
films could locally increase conductivity at the wall. Mean-field models considering 180° 
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domain walls in uniaxial ferroelectric semiconductors show that charging due to domain 
wall inclination significantly increases static conductivity [16J. Although energetically less 
favourable, such inclined domain walls can be present during polarisation switching, and 
can be locally stabilized within thin films by the presence of charged defects [T7j. We 
note that the dynamic nature of the domain wall conduction in BiFe03 was highlighted by 
recent measurements suggesting that microscopically irreversible distortions of the polari- 
sation structure are at its origin, and could be present in other ferroelectric domain walls 
[18] . Local probe studies in PZT have demonstrated a strong electroresistance effect, with 
conduction across thin films controlled by polarisation orientation [19] in agreement with 
theoretical predictions [20]. Essentially, the ferroelectric thin film with its adjoining elec- 
trodes (one of which can be the mobile metallic atomic force microscope (AFM) tip) is seen 
as two back-to-back Schottky diodes, with conduction dominated by interface or bulk lim- 
ited contributions depending on specific material parameters (dielectric response, mobility 
of charge carriers, presence of donor/trap states, metal and ferroelectric work functions and 
electrode geometry) [2TH23] . However, no observations of domain- wall-specific conduction 
have previously been reported in this material. 

All our measurements were carried out in ultra high vacuum on epitaxial PZT thin films 
grown on SrRu03 bottom electrodes on SrTiC>3 substrates. As described previously [2^1125] . 
the films show high crystalline and surface quality. The polarisation axis is perpendicular 
to the film plane, and all films are monodomain (up-polarised) as-grown. In these films 
we created 180° domain walls by switching the ferroelectric polarisation using a biased 
scanning AFM tip. The resulting domains were imaged by piezoresponse force microscopy 
(PFM) to determine the polarisation orientation, and by conductive-tip AFM (c-AFM) for 
local current maps. 

As shown in Fig. [T] for a 70 nm thick film, under a low negative DC tip bias we observe 
current signal at the domain walls, while the bulk of the sample remains insulating within our 
experimental resolution. In fact, as can be seen from the average signal profiles in Fig. [jje), 
the effective domain wall width measured by PFM (convoluted with the nominally ^45 nm 
radius tip) exactly corresponds to the lateral extent of this current signal. No correlation 
between the regions of high current in the domain walls and the morphology of the film 
was detected, and no features corresponding to the domain wall position were present in the 
surface topography. 
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FIG. 1: Domain wall conduction in PZT. Topography (a) showing an rms roughness of 0.4 nm. 
PFM phase (b) and amplitude (c) images of a square domain written with a positive tip bias, (d) 
c-AFM measurement at -1.5 V tip bias, showing an average current value of 25 pA at the position 
of domain walls, (e) Line profiles of c-AFM current and inverse PFM amplitude averaged over the 
left domain wall, as indicated by the dashed lines in (c) and (d). 

To better understand the microscopic mechanisms underlying the observed domain wall 
current (/), we investigated its dependence on the applied voltage (V). On stripe domains 
written in the same film, we carried out alternating PFM and c-AFM measurements with 
a stepwise increase of 0.125V for each iteration of the latter. As shown in Fig. [2](a), for 
negative tip voltages, we first detect current at -0.5 V. From the lack of any evolution of 
the ferroelectric domain structure imaged by PFM, it is clear that the observed current is 
not due simply to polarisation reversal and domain wall motion in the applied electric field 
(Fig. [2](b,c)). Moreover, static measurements at low voltages applied for 120 s likewise show 
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FIG. 2: Voltage dependence of domain wall conduction, (a) Current values from successive c- 
AFM scans with increasing tip bias, for sub-switching voltages from -0.50 V to -1.625 V, averaged 
over 120 scan lines, (b) -1.625 V c-AFM and (c) subsequent PFM measurement, (d) -2.25 V 
c-AFM measurement showing switching centers with current values on the order of 2-4 nA and 
(e) subsequent PFM measurement, demonstrating corresponding partial switching, (f) Average 
domain wall current as a function of the applied voltage with PFM measurements of polarisation 
hysteresis under applied voltage in inset. 

persistent domain wall current over 120 s with no change in the domain configuration (see 
Supplementary Materials). In fact, although displacement current associated with polari- 
sation switching is unambiguously seen at higher voltages (beyond -2.2 V), it is of much 
greater magnitude, and accompanied by changes in the domain configuration (Fig. [2](d,e)). 
Combining these data with similar measurements for positive tip bias, we can see that the 
I V characteristics of the domain wall are strongly non-linear and asymmetric (Fig. [2](f)). 
For positive bias, we do not observe domain wall current until much higher voltages (2.25 
V). We observe a similarly asymmetric polarisation switching hysteresis (inset) and relate 
these to the non-identical electrodes leading to different effective Schottky barrier heights. 

A more complete data set for the I-V characteristics was obtained from dynamic grid 
measurements on designated segments of the domain wall at different temperatures, as 
indicated in Fig. |3](a). For each 16 x 16 nm 2 sector of the grid, the voltage was swept 
from zero to V max and back to zero on the stationary tip, with increasing V max for each 
subsequent measurement on the same grid. As shown in Fig. [3](b), the I V characteristics 
averaged along the grid sectors containing the domain wall appear almost identical at 120 
and 150 K. At higher temperature, however, thermal activation is evident, with increasing 
current and decreasing conduction threshold voltage (Fig. |3](b) inset). More details of the 
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FIG. 3: Domain wall conduction at different temperatures, (a) PFM image of domain struc- 
ture with I-V measurement regions indicated by dashed lines, (b) Average current as a function 
of applied voltage, with the inset showing conduction thresholds (/ > 1.0 pA) for the different 
temperatures. 

measurements can be found in Supplementary Materials. 

We analyzed these data within the framework of different mechanisms includ- 
ing interface-limited conduction via Schottky thermionic emission (SE) log(//T 2 ) oc 
e/(k B T)(^eE/(4,7T£ e) - <j> B ) [26J or Fowler-Nordheim tunneling (FN) \og{I/E 2 ) oc 
—1/E [27] and bulk-limited conduction via Poole- Frenkel hopping (PF) \og(I/E) oc 



(—Utrap/ {kBT))(^J eE / (Attsqs)) or space charge limited conduction (SCL) / oc V n . In these 
expressions e is the electronic charge, T the temperature, E the electric field, k B the Boltz- 
mann constant, e the dielectric constant, £o the vacuum permittivity, (j) B the barrier height, 
and [/trap the trap energy. For SCL, different values of the exponent n are expected (n — 2 
for discrete traps, n > 2 with (n — 1) oc 1/T for traps distributed within the band gap) 
[32] . The electric field E across an insulating ferroelectric film is uniform for a parallel plate 
capacitor, but inhomogeneous in the case of an AFM tip, although a linear dependence 
on voltage can still be demonstrated (E oc V). In a semiconducting ferroelectrics parallel 
plate capacitor the presence of free charge, in particular screening polarisation near the film 
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surface, gives a maximum field E = y £ ^ c ^ (Y + Vu ~) m this region, with electronic 

transport limited primarily by the reverse-bias Schottky barrier at the film-electrode inter- 
face (abrupt junction approximation) [26], where N D is the dopant density, q the carrier 
charge, 6dc the static dielectric constant and Vu the resulting in-built bias. In the following 
discussion, we consider E oc V . Analysis of the data in the abrupt junction approximation 
yields unphysically high dopant concentration values (see Supplementary Materials). 

The high n values of the / oc V n fits of Fig. [3](a), and lack of obvious linearisation in 
n — 1 vs. 1/T suggest that SCL, either in the discrete trap limit or with a continuous trap 
distribution, is unlikely. However, representations of the data in normalised coordinates 
corresponding to FN, SE and PF all show reasonable linearisation (Fig. [i](b-d)). Very 
intense fields near the film surface, as modeled numerically in Supplementary Materials, 
would promote electron tunneling or emission into this region, suggesting an interface-limited 
mechanism in agreement with the significantly increased conduction when negative bias is 
applied to the tip rather than the bottom SrRu03 electrode. However, the obvious variation 
of the \og(I /V 2 ) vs. 1/V gradient with temperature (inset) suggests that FN, for which this 
quantity should be temperature-independent, does not adequately describe the data, while 
fits of log(7) vs. y/V in the SE regime yield unphysically low dielectric constant values, 
expected to be in the [e^ = 6]£dc — 400] range [21]. More realistic e values are obtained 
from the \og{I /V) vs. \/V fits in the PF regime, although these are still slightly smaller 
than the optical dielectric constant e^. Although our data preclude a quantitative Arrhenius 
analysis of the log(7) vs. 1/T dependence expected in this scenario, the almost identical I-V 
response observed at 120 and 150 K suggests more complex features. One possibility would 
be tunneling-assisted Poole- Frenkel conduction, previously observed in Hf02 thin films [28j, 
where tunneling/emission from the highly asymmetric electrodes gives rise to the diode-like 
I-V characteristics, but conduction within the domain wall itself proceeds by Poole-Frenkel 
hopping. 

To further explore the nature of the domain wall conduction, we also compared the current 
response for forward and reverse directions of voltage sweeps in the grid measurements 
carried out to progressively higher voltage V max: observing three distinct regimes. Initially, 
domain wall conduction with no hysteresis (calculated locally for each 16 x 16 nm 2 pixel 
as J /forward — Reverse) is detected (at -0.9 to -1.4 V in these measurements). Subsequently 
(-1.5 to -2.7 V), local hysteresis variations are observed (Fig. [5](c,d,g)) with no changes in 
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FIG. 4: Analysis of I—V characteristics for different conduction mechanisms, (a) linearisation in 
I oc V n (SCL), with (n - 1) vs. l/T in inset, (b) linearisation in I/V 2 vs. l/V (FN), with the 
gradient vs. T in inset, (c) linearisation I/T 2 vs. y/V (SE). (d) linearisation in I/V vs. W (PF). 

domain structure. These can be positive (higher current during forward voltage sweep) or 
negative (higher current during reverse voltage sweep) and vary randomly for consecutive 
voltage sweeps, including within the same grid sector, giving little to no hysteresis in the 
average I-V response, although negative hysteresis begins to dominate at the upper part of 
the range. Finally, upon reaching polarisation switching (beyond -2.8 V) we observe much 
higher current values, and strong negative hysteresis both at the domain walls and within 
the domains (Fig. [5](e,f,g, inset)), as expected for irreversible changes in the polarisation 
structure. 

Taken together, these data point to a thermally activated, domain-wall-specific phe- 
nomenon. As for conduction through thin ferroelectric films, asymmetric Schottky barriers 
and electrode tunneling into trap states near the interface appear to play a key role. However, 
there is clearly not enough trap density to support conduction across the ferroelectric film 
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FIG. 5: Cur rent- volt age hysteresis at domain walls. (a,c,e) C-AFM images of the domain wall 
conduction at V^aax = -1-4, -2.0, and -2.4 V. (b,d,f) Local hysteresis for voltage sweeps up to 
Vmax = -1-4, -2.0, and -2.4 V. (g) Average I-V domain wall response for the domain wall in the 
sub-switching regime: for the domain wall as a whole, no hysteresis is observed until voltages close 
to the switching threshold are reached. Inset shows the strong I-V hysteresis once polarisation at 
switching voltages. 

outside the domain walls. The latter must therefore act as a local trap reservoir, for which 
a range of activation energies could be envisaged. Such behaviour can be well correlated 
with the microscopic domain wall structure detected by TEM, with polarisation rotation 
near the film-electrode interface and small, apparently randomly spaced 1-2 unit cell steps 
throughout the length of the wall presenting maximally anti- aligned dipole moments fT5] . 
These charged domain wall segments necessitate screening, and would thus strongly increase 
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defect density at the wall. Defects such as oxygen vacancies could provide the necessary trap 
states for electronic tunneling from the electrode into the ferroelectric layer and conduction 
via subsequent hopping between trap states, and would be expected to show strong thermal 
activation. Rather than the uniform static conductivity predicted for inclined domain walls 
[T6] . we would expect significant local variations as a result of the microscopic structure - 
as are indeed observed - depending on both the density of charges/defects and the local 
connectivity of the charge regions around the domain wall steps. 

The lack of hysteresis at low bias suggests a strong pinning of the domain wall and accom- 
panying trap distribution in their initial metastable state, with predominantly conductive 
domain wall currents. At higher voltages, the microscopic structure of the domain wall could 
locally evolve at lengthscales smaller than the domain wall thickness/PFM resolution limit, 
with local hysteresis, different from the observations of global negative hysteresis through 
most of the conductive regime in BiFeOa [IS]- At high enough fields the depinning regime 
with large scale domain wall motion and reaccommodation of defects/screening charges ac- 
companying polarisation switching would result in a new metastable configuration. A crucial 
and surprising result in this scenario is the very high stability of the domain wall conduction 
over time. Although inclined 180° domain walls would thermodynamically be expected to 
relax to an elastically ideal straight configuration, leading to a decrease in the conductive 
signal, we observe persistent domain wall current over 120 s at low voltages, and the mea- 
sured current levels at a given domain wall segment remain stable when sampled over up to 4 
days at 295 K. We also note that in such PZT films, our previous measurements have always 
demonstrated very high stability of even 10-20 nm domains over weeks or months and at 
elevated temperatures [29], suggesting a high persistence of the configuration, presumably 
due to very effective defect-controlled pinning of the domain walls [30l EI] . This result is 
particularly encouraging for potential device applications, for which local conductive paths 
created using such domain walls should not deteriorate over time. 

In conclusion, our results clearly demonstrate conduction at 180° domain walls in ferro- 
electric tetragonal PZT, strongly suggesting that this phenomenon is far more general than 
has previously been shown. The scenario of defect segregation at partially charged domain 
walls could be present in a number of ferroelectric materials, and potentially controlled 
by strain and electrostatic boundary conditions, opening interesting possibilities for future 
research. In particular, studies looking for domain wall conduction in a large range of ferro- 
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electric materials should be envisaged to identify the best candidates for potential devices. 
The possibility of tuning the observed response by doping or oxygen vacancy engineering 
can also be explored. From a fundamental point of view, it would be interesting to compare 
the behaviour of domain walls in thin films and in single crystals, and to combine PFM, 
c-AFM and TEM studies, in particular to investigate the effects of intrinsic domain walls 
structure on the conduction. 

A. Methods 

60-70 nm PZT films on 30-40 nm metallic SrRu0 3 were epitaxially grown on single- 
crystal (001) SrTi03 substrates by off- axis radio- frequency magnetron sputtering. X-ray 
diffraction analysis showed relaxiation, with a c-axis lattice parameter values of ^4.16 A. 
AFM topography measurements showed 4 A rms surface roughness. 

PFM and c-AFM measurements were carried out in an Omicron FT beam deflection AFM 
under ultra-high vacuum (~3x 10 -10 mbar) using a Nanonis RC4 real-time controller and 
Bruker CoCr-coated MESP tips on samples introduced from ambient with no subsequent 
treatment. Typical scan size and rates were 3 /xm and 600 nm/s. For PFM, a 2 V, 1 kHz AC 
tip bias was applied, and PFM signal measured with an SRS SR830 lock-in amplifier (10 ms 
time constant). For c-AFM, a 0-10 V DC tip bias was applied, and current measured using 
an Omicron PRE E preamplifier with a ±500 fA noise level. Conduction threshold values 
depend on tip condition, with ^400 mV for the best tips. 
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